MUiJi'i:iiiii\,  Ojtiijii'v^i-oWlA  ybSJ-ib-SOOS 


A 


DEPARTMENT  OF  OCEAN  ENGINEERING 

MASSACHUSEnS  INSTITUTE  OF  TECHNOLOGY 
CAMBRIDGE.  MASSACHUSEHS     02139 


;,f%| 


COMPARISON  OF  ENTROPY  GENERATION  AND 
CONVENTIONAL  DESIGN  METHODS  FOR  HEAT  EXCHANGERS 

by 
David  Samuel  Herbein 


XIII-A 


June    1987 


j^^--v<g>«^.. 


■Si' ■■■■7  ;3Tt'^"  -  ••>       -^s^. 


^'^-;--: 


(65*^^ 'jK 


-#^';^^V; 


^fe^^ 


.-wf--'*sa=^     t-riJ 


-^1^,  -^1 


..'"""^^  -^^ 


COMPARISON  OF  ENTROPY  GENERATION  AND 
CONVENTIONAL  DESIGN  METHODS  FOR  HEAT  EXCHANGERS 

by 

David  Samuel  Herbein 

B.S.  Marine  Engineering 

United  States  Naval  Academy 

(1975) 

Submitted  to  the  Department  of  Ocean  Engineering 

in  Partial  Fulfillment  of  the  Requirements  for  the 

Degrees  of 

NAVAL  ENGINEER 

and 

MASTER  OF  SCIENCE  IN  NAVAL  ARCHITECTURE  AND 

MARINE  ENGINEERING 

at  the 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 

June,  1987 
©David  Samuel  Herbein,  1987 

The  author  hereby  grants  to  M.I.T.  and  to  the  U.S.  Government  permission 
to  reproduce  and  to  distribute  copies  of  this  thesis  document  in  whole  or  in 
part.  / 


COMPARISON  OF  ENTROPY  GENERATION  AND  CONVENTIONAL 
DESIGN  METHODS  FOR  HEAT  EXCHANGERS 

by 
David  Samuel  Herbein 

Submitted  to  the  Ocean  Engineering  Department  on  May  8,  1987  in  partial 
fulfillment  of  the  requirements  for  the  degrees  of  Naval  Engineer  and  Master 
of  Science  in  Naval  Architecture  and  Marine  Engineering. 

Abstract 

The  design  of  heat  exchangers  traditionally  focuses  on  the  known  con- 
straints of  the  problem  such  as  inlet  and  outlet  temperatures,  flow  rates,  and 
pressure  drops.  This  leads  mainly  to  a  sizing  problem  where  the  designer 
must  select  surfaces,  flow  configuration,  and  materials  to  meet  the  mini- 
mum design  objectives.  An  alternate  approach  based  on  an  acceptable  level 
of  thermodynamic  irreversibility  (entropy  generation)  has  been  proposed. 
When  the  entropy  generation  level  has  been  set,  the  geometric  parameters 
of  the  heat  exchanger  can  be  determined.  The  design  of  a  plate-fin  type, 
gas-to-gas  recuperator  for  a  regenerative  open  Brayton  cycle  has  been  used 
as  a  demonstrative  device.  The  resulting  heat  exchanger  designs  are  then 
examined  to  determine  what  caused  the  differences  and  why  either  method 
should  be  preferred  over  the  other. 
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Chapter  1 

INTRODUCTION 

1.1  Background 

Energy  conservation  is  a  topic  that  has  recently  received  considerable  atten- 
tion. Available  work  has  been  recognized  as  a  valuable  commodity  and  its 
destruction  in  seemingly  otherwise  efficient  engineering  processes  has  also 
seen  a  renewed  increase  in  interest.  The  use  of  second  law  analysis  and  ther- 
modynamic irreversibility  minimization  have  been  proposed  as  techniques 
that  should  become  an  integral  part  of  the  design  of  engineering  processes 
and  components  [l].  It  is  not  apparent  what  the  relationship  of  traditional 
methods  of  design  and  these  irreversibility  minimization  techniques  will  be. 

1.2  Conventional  Design  Procedures 

The  design  of  a  component,  such  as  a  heat  exchanger,  generally  involves  the 
sizing  of  that  component  to  meet  specified  performance  parameters  within 
known  constraints.  The  heat  exchanger  area  and  volume  are  usually  de- 
signed to  be  the  minimum  required  to  meet  the  specifications  as  this  will 
also  usually  be  the  best  design  from  an  ecomonic  aspect.  The  fluid  inlet  and 
outlet  temperatures,  flow  rates,  and  pressure  drops  are  usually  specified  and 
it  is  the  designer's  task  to  determine  construction  type,  flow  arrangement, 
materials  and  surfaces  to  meet  those  requirements  [2].  The  effect  of  the 
component  design  on  the  overall  system  performance  should  also  be  consid- 
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ered  to  ensure  the  tradeoffs  performed  within  its  design  are  still  valid  when 
integrated  into  the  system. 

1.3      Irreversibility  Minimization 

Any  heat  transfer  process  is  generally  accompanied  by  thermodynamic  ir- 
reversibility or  entropy  generation.  The  entropy  generated  is  in  direct  pro- 
portion to  the  amount  of  useful  work  dissipated  in  that  process.  If  this 
irreversibility  can  be  minimized  within  a  particular  component  of  a  power 
cycle,  the  useful  power  output  of  the  cycle  should  increase  [3].  Bejan  has 
shown  how  the  entropy  generation  rate  can  be  reduced  in  a  counterflow  gas- 
to-gas  heat  exchanger  [4],  and  it  is  this  method  that  will  be  used  to  create 
designs  that  will  be  compared  to  conventionally  designed  heat  exchangers, 
and  their  impact  on  overall  cycle  performance. 
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Chapter  2 

SELECTION  OF 
BRAYTON  CYCLE 
PARAMETERS 

2.1      Cycle  Components 

A  regenerative  open  Brayton  cycle  was  chosen  as  the  vehicle  to  test  the 
different  design  methods  for  the  recuperator.  A  schematic  of  the  cycle  com- 
ponents is  shown  in  Figure  2.1. 


Turbine 


Compressor 


Regenerative    Cycle    Schematic 
Figure    2.1 
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A  compressor  polytropic  efficiency  was  assumed  and  used  throughout 
all  calculations,  rjcp  =  0.95.  Similarly,  a  turbine  polytropic  efficiency  was 
assumed,  rjct  =  0.90.  For  ease  of  calculations,  only  one  fluid  was  used, 
i.e.,  air  was  modeled  as  a  perfect  gas  with  Cp  =  1.0  kj/kg  —°K  and  R  = 
0.287  kj/kg— °  K  both  assumed  constant.  The  combustor,  then,  was  treated 
as  a  perfect  heat  transfer  device,  but  was  not  an  injection  point  for  fuel. 

2.2      Cycle  Parameters 

A  reasonable  compressor  inlet  temperature  was  selected,  Ti  =  300° K,  along 
with  an  inlet  pressure  of  Pi  =  1.126  x  lO^iV/m^.  A  nominal  value  of  turbine 
inlet  temperature  was  also  selected,  T4  =  1300°  ii".  The  selection  of  an 
appropriate  pressure  ratio  required  the  following  consideration. 

The  addition  of  a  recuperator  to  a  simple  open  Brayton  cycle  causes  a 
shift  in  the  pressure  ratio  for  maximum  cycle  efficiency  from  that  of  the  non- 
regenerative  cycle.  A  heat  exchanger  effectiveness  of  e  =  0.8  was  selected 
and  cycle  efficiency  without  losses  was  calculated  for  pressure  ratios  varying 
from  4  to  10.  A  sample  calculation  is  contained  in  Appendix  A,  results  are 
plotted  in  Figure  2.2. 
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Based  on  these  results,  a  pressure  ratio  of  P^jPx  =  5  was  selected.  With 
this  pressure  ratio  determined,  the  maximum  cycle  efficiency  without  losses, 
r7oo,  could  be  calculated.  This  quantity  will  be  used  for  comparison  purposes 
and  was  computed  as  r]^  =  0.5760. 

2.3      Heat  Exchanger  Characteristics 

A  plate-fin  heat  exchanger  surface  was  selected  from  those  cataloged  in  ref- 
erence [5].  Its  characteristics  and  dimensions  are  shown  in  Figure  2.3.  It  has 
an  uncomplicated  geometry  and  a  reasonably  wide  range  of  Reynolds  num- 
bers in  which  consistent  results  can  be  expected.  The  material  is  aluminum, 
with  a  thermal  conductivity  o(kf  =  190  w/m  —°  K. 


Fin  pitch  =  2.0  per  in  =  78.74  per  m 

Plate  spacing,  6  =  0.750  in  =  19.05  x  10"' m 

Fin  length  =  12.0  in  =  304.8  x  10-'m 

Flow  passage  hydraulic  diameter,  Ar^  =  0.0474  ft  =  14.453  x  10'^m 

Fin  metal  thickness  =  0.032  in,  aluminum  =  0.813  x  10'^m 

Total  heat  transfer  area/volume  between  plates,  0  =  76.1  ft'/ft-'  =  249.672  m'/m-" 

Fin  area/total  area  =  0.606 


Plain   Plate-Fin   Sur-^ace   3.0 
Figure  ^3 
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Chapter  3 

CONVENTIONAL 
DESIGN  METHOD 

3.1  Introduction 

This  chapter  outlines  the  conventional  design  procedure  that  was  used  to 
size  the  single-pass,  counterflow  gas-to-gas  heat  exchanger.  This  procedure 
was  performed  for  various  combinations  of  heat  exchanger  effectiveness  and 
specified  pressure  drops.  A  complete  numerical  example  is  provided  in  Ap- 
pendix B.  Final  results  for  this  conventional  method  are  shown  in  Table 
3.1. 

3.2  Conventional  Methodology 

For  a  single-pass  counterflow  design,  when  the  core  dimensions  on  one  side 
are  fixed,  the  dimensions  on  the  other  side  are  also  known.  This  means  the 
design  is  driven  by  the  side  that  has  the  more  stringent  AP/P  requirement. 
The  method  of  determining  the  controlling  side  is  contained  in  Appendix 
B.  In  every  case  considered,  the  hot  side  of  the  heat  exchanger  was  the 
controlling  side  and  so  its  requirements  drove  the  design. 

The  following  is  a  step-by-step  heat  exchanger  sizing  design  procedure 
which  closely  follows  that  outlined  in  reference  [2]. 

1.  For  the  given  heat  exchanger  effectiveness,  determine  the  fluid  outlet 
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temperatures.  Calculate  the  fluid  mean  temperature  on  each  side  and 
evaluate  fluid  physical  properties  /),-,  po,  (l/p)m,  and  fx. 

2.  Determine  Ntu  for  the  exchanger  and  then  Ntu  for  each  side.  The  in- 
fluence of  longitudinal  heat  conduction  is  ignored  in  this  first  iteration 
of  design. 

3.  Estimate  hot  and  cold  side  pressure  drops,  select  an  appropriate  Njie^ 
and  then  a  value  of  j/f  from  Figure  2.3. 

4.  Calculate  mass  velocity,  G,  from  information  in  steps  1-3  and  the  cor- 
responding value  of  AP/P. 

5.  Calculate  Nrc  and  determine  values  of  j  and  f  from  Figure  2.3. 

6.  Compute  heat  transfer  coefficient,  h;  temperature  effectiveness  of  the 
fins,  r/y;  and  the  total  surface  temperature  effiectiveness,  rjo. 

7.  Calculate  heat  transfer  area,  A;  minimum  free  flow  area,  Ao]  heat  ex- 
changer frontal  area,  A/ri  flow  length,  L;  and  heat  exchanger  volume, 
Vol. 

8.  Compute  AP/P  from  known  conditions  and  calculated  parameters. 

3.3      Analysis 

For  each  case  show  in  Table  3.1,  cycle  efficiency  with  losses,  77,  was  calculated 
and  the  cycle  efficiency  degrade,  rjoo  —  »7j  is  listed  for  each  case.  A  correlation 
between  heat  exchanger  effectiveness,  e,  total  allowed  pressure  drop,  k,  cycle 
efficiency  degrade,  and  heat  exchanger  volume  was  sought. 

The  relationship  between  cycle  efficiency  degrade  and  effectiveness  at 
various  levels  of  k  is  shown  in  Figure  3.1.  The  increased  level  of  degrade 
at  lower  values  of  effectiveness  is  as  expected,  and  for  a  particular  value  of 
eff'ectiveness,  the  degrade  increases  with  increasing  pressure  drops. 

The  effects  of  varying  effectiveness  and  pressure  drops  on  heat  exchanger 
volume  is  depicted  in  Figure  3.2.    The  higher  the  effectiveness,  the  longer 
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and  therefore,  larger  the  heat  exchanger.    Similarly,  the  greater  pressure 
drop  for  a  particular  effectiveness  also  increases  the  volume. 

Figure  3.3  was  derived  by  extracting  values  of  effectiveness  at  nominal 
values  of  r^oo  —  ^  and  plotting  them  on  the  appropriate  k  curve  of  Figure  3.2. 
This  resulted  in  curves  of  constant  r/oo  —  H  that  have  a  minimum  volume 
at  a  particular  value  of  k  and  e.  It  is  emphasized  that  these  "minimum" 
volumes  were  obtained  at  constant  values  of  cycle  degrade  with  varying 
pressure  drops  and  effectiveness.  These  "minimum"  volumes  are  plotted 
against  effectiveness  in  Figure  3.4  with  the  respective  k  values  annotated  on 
the  graph. 
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Case  # 

e 

k 

noo  -  n 

L(m) 

Vol  (m^) 

C-1 

1-0 

0.02 

0.0043 

— 

— 

C-2 

1.0 

0.04 

0.0089 

— 

— 

C-3 

1.0 

0.06 

0.0137 

— 

— 

C-4 

1.0 

0.08 

0.0188 

— 

— 

C-5 

0.9 

0.0 

0.0449 

— 

— 

C-6 

0.9 

0.0240 

0.0498 

13.69 

3.69 

C-7 

0.9 

0.0388 

0.0550 

13.82 

2.64 

C-8 

0.9 

0.0558 

0.0603 

14.04 

2.19 

C-9 

0.9 

0.0739 

0.0660 

14.65 

1.98 

C-10 

0.8 

0.0 

0.0833 

— 

— 

C-11 

0.8 

0.0192 

0.0886 

6.15 

1.12 

C-12 

0.8 

0.0367 

0.0941 

6.56 

0.85 

C-13 

0.8 

0.0550 

0.0999 

6.72 

0.71 

C-14 

0.8 

0.0740 

0.1060 

7.15 

0.65 

C-15 

0.7 

0.0 

0.1165 

— 

— 

C-16 

0.7 

0.0188 

0.1221 

3.75 

0.53 

C-17 

0.7 

0.0365 

0.1278 

3.95 

0.39 

C-18 

0.7 

0.0565 

0.1338 

4.17 

0.34 

C-19 

0.7 

0.0734 

0.1401 

4.32 

0.31 

C-20 

0.6 

0.0 

0.1455 

— 

— 

C-21 

0.6 

0.0186 

0.1512 

2.52 

0.29 

C-22 

0.6 

0.0376 

0.1571 

2.69 

0.22 

C-23 

0.6 

0.0557 

0.1632 

2.85 

0.19 

C-24 

0.6 

0.0764 

0.1695 

2.98 

0.17 

Table  3.1:  Conventional  Method  Results 
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Chapter  4 

MINIMUM  ENTROPY 
DESIGN  METHOD 

4.1  Introduction 

The  use  of  irreversibility  concepts  in  the  design  of  the  gas-to-gas  counterflow 
heat  exchanger  will  be  presented  in  this  chapter.  For  the  same  mass  veloci- 
ties that  were  used  in  Chapter  3,  the  heat  exchanger  geometric  parameters 
can  be  determined  to  arrive  at  minimum  levels  of  entropy  generation.  The 
outputs  of  this  design  procedure  are  heat  exchanger  effectiveness  and  the 
pressure  loss  suffered  on  each  side.  For  a  more  detailed  explanation  of  this 
procedure,  the  reader  is  directed  to  references  [4]  and  [6]. 

4.2  Background 

Since  the  method  deals  with  the  entropy  generation  within  the  recuperator, 
it  is  appropriate  to  begin  with  an  entropy  flux  analysis  on  the  control  volume 
of  the  recuperator  in  Figure  4.1.  The  heat  transfer  from  the  outer  walls  is 
assumed  to  be  negligible.  The  entropy  generation  rate  can  then  be  written 
as 

S  =  Thi{Sout  -  Sin)l  +  rh2{Sout  -  Sm)2  (4.1) 
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Strean   1 


/\  \/ 


Stream   2 


T, 


Recuperator"  Control    Volume 
Figure   4.1 


Substituting  expressions  for  the  ideal  gas  entropy  changes  based  on  pressure 
and  temperature  relationships,  this  becomes 


S  =  CminMTiout/Ti)  +  {R/Cp)iln{Pi/Piout)] 
+CmaxMT2out/T2)  +  {R/Cp)2ln{P2/ Piout)] 

The  number  of  entropy  production  units,  A^^,  is  defined  as 


(4.2) 


(4.3) 


Writing  expressions  for  the  first  law  of  thermodynamics  of  the  control  vol- 
ume, 

Cmin{Tl  -  Tiout)  +  Crmn{T2  "  T2ont)  =  0  (4.4) 

and  heat  exchanger  effectiveness. 


€={Tiout-Ti)/{T2-Ti) 


(4.5) 
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permits  tbe  elimination  of  Tiout  and  T2out  in  equation  (4.2)  which  allows  the 
rate  of  entropy  production  formula  to  be  nondimensionalized  as 


AT  ^min  f 

N.  =  — In 


X.e||-1 


+  ln 


1  - 


7\ 
T2 


M  (1  -  ^)  (^.6) 

'P/l  ^  ^      /  1         \^pj  2         ^  ^    ^2 

A  special  form  of  this  last  expression  can  be  obtained  in  the  case  of  nearly 
ideal  heat  exchangers.  When  the  stream-to-stream  AT's  and  frictional  AP's 
are  held  small,  two  inequalities  can  be  applied. 


1-e     << 1 


(4.7) 


(AP/P)i,2      «  1 
Applying  these  and  expressing  e  in  terms  of  Ntu  yields 


Ns     =^ ln—  +  in 

^max       -1 1 


1  - 


1- 


T2 


(4.8) 


+ 


1  - 


^min   1 

(>  -  ^)' 

exp 

'-N 

^u(l- 

^max  j  , 

^max  J 

1  - 

Cmax     \               T2  J 

1    - 

^max 

exp 

-iVtu 

f  1  -  ^^^ 

(4.9) 


The  first  two  terms  can  be  characterized  as  a  contribution  from  capacity 
rate  imbalance,  the  third  represents  contribution  due  to  finite  Ntu,  and  the 
last  two  show  fluid  friction  eff'ects. 

Requiring  the  streams  to  be  balanced  and  applying  the  calculus  of  limits 
as  Cmin  ~*  Omax,  equation  (4.9)  can  be  reduced  to 


Ns  ci In—     +     In 


1  - 


^mir 


C„ 


1-^ 

T2 
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When  the  overall  Nm  is  expressed  in  terms  of  Ntu  for  each  side  of  the  heat 
transfer  surface, 

^      __        ^         ,     ^min       1  f4  1  ll 

Ntu         Ntul         Cmax  Ntu2 

the  number  of  entropy  generation  units  can  be  divided  into  three  contribu- 
tions 

NsC:^  Ns  imbalance  +  Nsl  +  Ns2  (4.12) 

Applying  the  limit  Cmin  ~*  Cmax,  the  imbalance  terms  of  equation  (4.10) 
can  be  written  as 

Nsimbalance  =      7; M   I  ;p  -  1  -  7; In—  (4.13) 

which  will  clearly  vanish  when  Cmin  =  Cmax-  Since  this  is  the  case  for  the 
design  here,  the  number  of  entropy  production  units  for  each  side  can  be 
determined  from 

4.3      Design  Considerations 

The  number  of  entropy  generation  units  that  have  been  deemed  acceptable 
for  a  particular  design  is  a  function  of  the  hydrauhc  radius,  minimum  free 
flow  area  and  flow  length.  It  is  necessary,  then,  to  express  Ng  in  terms  of 
these  and  other  known  flow  parameters.  Using  the  definitions  of  Ntu.  and 
friction  factor  for  each  side 

iVtai,2    =     {L/rh)Nst  (4.16) 

(AP/P)i.2     =     fiL/rh)Gy{2pP)  (4.17) 
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in  equations  (4.14)  and  (4.15)  leads  to 


where 


Nsi,2  =  tttTTaT  +  ^Bfi^/'^)3'  (4.18) 


ai  =  1 

a2  =  Cmin/Cmax  (4.19) 

"1  ^^  ^min/^max 

62  =  1  (4.20) 

Si  -  (i?/Cp)i 

B2  =  {R/C,)2  (4.21) 


'T2  Ti 

r      = 


2 


(4.22) 


g      =      G/{2pP)2  (4.23) 

As  can  be  seen  in  equation  (4.18),  the  ratio  L/rh  performs  a  trade-off  func- 
tion, i.e.,  for  a  fixed  g  and  Nrc,  there  will  be  an  optimum  L/rh  which  results 
in  a  minimum  iVsi_2.  When  the  function  Nsi^2  is  minimized,  the  {L/rh)opt 
is  given  by 

4.4      Minimum  Entropy  Design  Methodology 

The  method  used  is  based  on  optimizing  {L/rh).  Since  rh  is  set  by  the  selec- 
tion of  the  plate-fin  surface,  this  yields  a  length  that  will  produce  minimum 
entropy  for  the  given  hydraulic  radius. 
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The  following  outlines  the  methodology  employed.  A  complete  numerical 
example  is  shown  in  Appendix  C.  Final  results  for  all  cases  are  tabulated  in 
Table  4.1. 

1.  For  the  cases  listed  in  Table  3.1  that  resulted  in  a  heat  exchanger 
volume,  use  calculated  G  and  Njie  to  determine  {L/rh)opt- 

2.  Compute  resultant  heat  exchanger  effectiveness  and  hot  and  cold  side 
pressure  drops. 

3.  Calculate  N..,  N.,  and  N.,  ,  ,. 

"h  >  'c  '  "total 

4.  Calculate  the  heat  exchanger  volume  for  each  case. 

5.  Determine  the  cycle  efficiency  degrade,  r/oo  —  V- 

4.5      Analysis 

Figure  4.2  shows  the  nearly  linear  relationship  of  cycle  efficiency  degrade  and 
heat  exchanger  effectiveness.  There  is  not  the  clear  definition  of  different 
values  of  pressure  drop  as  wa.s  observed  for  the  conventional  method  in 
Figure  3.1.  In  fact,  there  is  even  a  wider  range  of  A:  values  here,  from 
k  =  0.056  for  case  E-6  to  A:  =  0.269  for  case  E-24,  and  yet  the  results  are 
almost  linear  across  that  range. 

The  volume  plotted  in  Figure  4.3  is  that  volume  obtained  by  minimizing 
entropy  generation  in  the  heat  exchanger  and  should  not  be  interpreted  as 
a  minimum  volume,  but  rather  as  an  "optimum"  volume  for  minimizing 
entropy  generation. 

To  achieve  a  better  ba^is  for  comparison  of  the  two  methods,  the  entropy 
generated  by  the  various  designs  must  be  discussed. 
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Case  # 

Lopt{m) 

Vol(m2) 

k 

€ 

'7oo  -  ^ 

Ns 

E-6 

32.03 

8.63 

0.0563 

0.955 

0.0347 

0.0325 

E-7 

25.86 

4.93 

0.0727 

0.945 

0.0440 

0.0417 

E-8 

22.04 

3.44      ■ 

0.0876 

0.935 

0.0528 

0.0503 

E-9 

20.12 

2.72 

0.1014 

0.926 

0.0606 

0.0583 

E-11 

24.35 

4.46 

0.07658 

0.942 

0.0461 

0.0441 

E-12 

19.08 

2.46 

0.1069 

0.922 

0.0641 

0.0614 

E-13 

16.11 

1.07 

0.1320 

0.908 

0.0784 

0.0760 

E-14 

15.08 

1.38 

0.1559 

0.896 

0.0917 

0.0897 

E-16 

19.80 

2.77 

0.996 

0.927 

0.0596 

0.0573 

E-17 

15.40 

1.53 

0.1422 

0.902 

0.0844 

0.0823 

E-18 

13.19 

1.07 

0.1778 

0.884 

0.105 

0.1028 

E-19 

12.11 

0.86 

0.2058 

0.872 

0.1211 

0.1189 

E-21 

16.98 

1.93 

0.1252 

0.913 

0.0741 

0.0718 

E-22 

12.87 

1.04 

0.1799 

0.882 

0.1067 

0.1048 

E-23 

11.52 

0.76 

0.2251 

0.862 

0.1335 

0.1309 

E-24 

10.47 

0.60 

0.2685 

0.847 

0.1596 

0.1541 

Table  4.1:  Entropy  Generation  Method  Results 
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Chapter  5 

COMPARISONS  AND 
RESULTS 

5.1  Introduction 

To  this  point,  the  minimum  entropy  method  results  have  been  compared  to 
the  conventional  method  results  in  terms  of  conventional  method  standards, 
i.e.,  pressure  drops,  efficiencies,  and  effectiveness.  It  is  necessary  to  compare 
these  methods  in  terms  of  the  minimum  entropy  method,  i.e.,  the  entropy 
generated  in  each  design. 

5.2  Conventional  Method  Entropy  Generation 

The  number  of  entropy  generation  units  for  each  design  can  be  computed 
from  previously  determined  information  and  equations  (4.14)  and  (4.15). 
An  example  of  these  calculations  is  contained  in  Appendix  D.  Results  for 
the  conventional  causes  are  listed  in  Table  5.1.  The  large  deviation  in  N3  and 
nc»  —  V  for  cases  C-21  through  C-24  are  a  result  of  violating  the  assumption 
that  1  —  €  <<  1.  These  conventional  cases  will  be  discarded  from  further 
consideration.  However,  the  minimum  entropy  cases  at  corresponding  mass 
velocities  are  still  valid. 

When  the  entropy  generation  units  are  plotted  against  the  degrade  in 
cycle  efficiency,  Figure  5.1,  the  dependence  on  pressure  drop  is  still  present. 
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The  grouping  of  data  points  is  caused  by  the  discrete  changes  in  heat  ex- 
changer effectiveness,  decreasing  from  left  to  right. 

Figure  5.2  was  created  by  using  the  same  values  of  rjoo  —  rj  that  were  used 
for  Figure  3.4.  The  value  of  N3  at  each  of  these  cycle  efficiency  degrades 
was  read  off  at  a  particular  value  of  k.  Plotting  these  "minimum"  volumes 
against  entropy  generation  rate  exhibits  the  general  trend  that  the  larger 
the  volume  of  the  heat  exchanger,  the  lower  the  total  entropy  generated. 
This  is  in  consonance  with  the  results  predicted  by  reference  [4]. 

5.3  Minimum  Entropy  Generation 

The  cycle  efficiency  degrade  vs.  entropy  generated  for  the  minimum  entropy 
method  is  similar  in  overall  shape  to  that  of  the  conventional  method,  but, 
as  Figure  5.3  shows,  the  relationship  is  nearly  linear  and  does  not  show  the 
discrete  pressure  drop  differences. 

Figure  5.4  shows  that  the  "optimum"  volumes  determined  hy{L/rh)opt 
follow  essentially  the  same  path  as  that  of  the  conventional  method.  For  a 
direct  comparison,  Figures  5.2  and  5.4  are  plotted  together  on  Figure  5.5. 

To  ensure  that  the  {L/rh)opt  expression,  (4.24),  was  giving  minimum 
entropy  results,  a  series  of  calculations  was  performed  at  [L/rh]  values  on 
both  sides  of  the  optimum  for  case  E-9.  Figure  5.6  shows  that  the  minimum 
entropy  solution  was,  in  fact,  being  determined.  It  is  also  noted  that  the 
magnitudes  of  N3  and  r^oo  —  »?  are  essentially  the  same  over  this  range  of 
heat  exchanger  effectiveness. 

5.4  Comparisons 

The  calculations  performed  on  case  E-9  to  prove  minimization  also  yielded 
heat  exchanger  volumes  at  a  constant  mass  velocity  G.  These  additional 
data  points  are  plotted  on  Figure  5.7.  The  dash-dot  Une  connecting  them 
is  a  locus  of  recuperator  volumes  along  a  constant  G  path.  This  shows  why 
the  volumes  at  minimum  A^^  or  r/oo  —  r}  along  a  constant  G  path  are  larger 
that  the  miminum  volume  of  a  given  Ns  or  7700  —  r). 
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Figure  5.8  shows  a  comparison  between  conventional  "minimum"  vol- 
umes and  minimum  entropy  "optimum"  volumes.  The  minimum  entropy 
method  yields  volumes  that  are  larger  than  the  conventional  method  for  the 
same  cycle  efficiency  degrade. 
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Case  # 

^co  -  ^ 

Ns 

C-1 

0.0043 

— 

C-2 

0.0089 

— 

C-3 

0.0137 

— 

C-4 

0.0188 

— 

C-5 

0.0449 

— 

C-6 

0.0498 

0.0445 

C-7 

0.0550 

0.0494 

C-8 

0.0603 

0.0550 

C-9 

0.0660 

0.0610 

C-10 

0.0833 

— 

C-11 

0.0886 

0.0900 

C-12 

0.0941 

0.0967 

C-13 

0.0999 

0.1036 

C-14 

0.1060 

0.1108 

C-15 

0.1165 

— 

C-16 

0.1221 

0.1505 

C-17 

0.1278 

0.1584 

C-18 

0.1338 

0.1669 

C-19 

0.1401 

0.1748 

C-20 

0.1455 

— 

C-21 

0.1512 

0.2308 

C-22 

0.1571 

0.2405 

C-23 

0.1632 

0.2501 

C-24 

0.1695 

0.2607 

Table  5.1:  Conventional  Method  Entropy  Generation 
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Chapter  6 

CONCLUSIONS 


6.1  Discussion 

The  nature  of  the  results  and  conclusions  drawn  for  the  example  cycle  would 
be  the  same  for  other  cycles  and  recuperator  surfaces.  For  high  effective- 
ness recuperators,  the  counterflow  arrangement  is  dictated.  Here  with  the 
same  finned  surface  on  both  sides,  the  (AP/P)c/(AP/P)/i  ~  [Ph/PcY  or 
approximately  1/25.  If  a  closer  plate  spacing  were  used  on  the  clean  air  cold 
side  more  AP/P  would  appear  on  the  cold  side  for  any  total  HP / P.  This 
would  result  in  even  smaller  minimum  volumes  at  any  given  Ns  or  r^oo  -  '?• 

6.2  Conclusions 

In  selecting  a  gas  turbine  recuperator,  essentially  the  same  results  are  ob- 
tained either  by  the  conventional  method  of  minimizing  degradation  of  cycle 
efficiency,  rjoo  —  r^,  or  by  minimizing  entropy  generation,  Ng. 

The  selection  of  an  "optimum"  recuperator  requires  that  care  be  exer- 
cised in  the  selection  of  the  path  along  which  the  optimization  is  calculated. 
It  has  been  shown  that  the  volume  of  the  recuperator  for  minimum  Ng  or 
»7oo  ~  ^  along  a  constant  mass  velocity  path  is  much  greater  than  the  true 
minimum  volume  at  any  r]^  —  r)  or  N3. 
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Appendix  A 

Sample  Calculation  for 
Effect  of  Recuperator  on 
Cycle  Efficiency 


e  =  O.S,r=- 


»7oo  = 


<Vnet  tVt  +  iUc 


Q  Q 

Cjt  =  m(/i4  -  h^)  =  mCp{T4  -  T5) 
Cjc  =  m{hi  -  /12)  =  rhCp{Ti  -  T2) 


Turbine: 


Q  =  m{h4  -  /13)  =  Thcp{T4  -  T3) 


T4    _    mN^^'-'^^V) 
T5     ~     \Ps. 


n  = 


-1\T-1 
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Ts     = 


■(6)(-^)('^)' 


1300 


=  820.1°fc 


Compressor: 

T2 

(P2\  Inlc)^) 

r,   =   Ti    -f 


P2\  (nlc)h) 


Recuperator: 


T2   =    {300)  {e)TmrT^  =  514.2°  k 


e     = 


T3-T2 
T5-T2 


Tz    =     {e){Ts-T2)  +  T2 


Ts     =     (.8)(820.1  -  514.2)  +  514.2  =  695.2°A: 

Cycle  Efficiency: 

(T4  -  Ts)  +  (Ti  -  T2)       (1300  -  820.1)  +  (300-514.2) 
T]qo  =  = = =  -r-rr =  U.4D0i 


T4-T3 


1300  -  695.2 
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Appendix  B 

Conventional  Design 
Method  Numerical  Example 


Before  calculations  can  be  performed,  two  points  need  to  be  clarified:  1)  the 
method  of  determining  how  loss  of  turbine  work  due  to  pressure  drops  in 
the  system  is  distributed  and;  2)  the  method  for  determining  which  side  of 
the  heat  exchanger  has  the  more  stringent  AP/P  requirement. 
Turbine  work  is  given  by  the  expression 


Wt 


rnr)tCpTh 


1- 


•*  c 


Wt 


=      1- 


mr}tCpTh 
Expanding  the  last  term, 


p.  +  APcY  fPcV 


APc 


Wt 


rnVt^pTh 


-  '-^t 


PhJ 


1  +  ^ 

■*    C 

1_   AP^ 


Ph 


Wt 


^       1-       ^ 


V  Ph 


l  +  7l  ^7^+-^)  + 


rhrjtCpTh  \Ph 

so  that  loss  of  turbine  work  can  be  expressed  as 
Ati; 


rnVtCpTh 


._,,^Vf^^^^c 


PhJ    V  Ph 
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and  it  is  clear  from  Figure  B.l  that  S^-  determines  the  loss  of  turbine 
work  irregardless  of  how  it  is  distributed.  Introducing  the  definitions  r  =  -W^ 


ID 
< 
UJ 

•z. 
u 


specific  entropy 

Temperature    -    Entropy    Diagran 
Figure    B.l 


and  Q  =  \^\       ,  it  can  be  shown  that  for  a  counterflow  exchanger  where 


Ac  =  Ah  and  Lc  =  Lh, 


APh/Pk        9 


APc/Pc        r2 
Using  a  slightly  different  form  of  the  last  turbine  work  expression, 


Awt 


rnVtCpTh 


1^ 


+  1 


AP. 


Awt 


rnrjtCpTh 


1    f9\AP, 
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Awt     _     _     1  {9/r^  +  1)  APc 


wt  'ri{l-l/ri)    P, 


Ati^     ^     _    {e/r^  +  1)  AP, 
wt  '^     ri  -1       P, 

Introducing  -p-^  +  -p-^  =  k,  expressions  for  hot  and  cold  side  pressure 
drops  are, 

APf,  k{9/r^) 


Ph  1  +  0/r^ 

APg       _  k 

P,-      ~     l  +  ^/r2 

The  numbered  steps  that  follow  correspond  to  those  outlined  in  Section 
3.2.  The  input  design  parameters  used  for  the  example  calculations  were 
e  =  0.8  and  k  =  0.04. 

B.l      Temperatures  and  Fluid  Properties 


r,     =     300°  A: 


1^2  /  ^2  \  ''pcT' 


SO  that 


T,    =    Txl'^^^ 


T2     =     300(5)  T^^^TTT^  =  486.8°  fc 


^     =     ^(1  -  A:)  =  5(1  -  .04)  =  4. 
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so  that 


T.    = 


P4yp'(V-) 
pj 


(a) 


npt{hr) 


T4 


-1 


Ts     = 


£      = 


(4.8)(-^)(^) 


1300 


-1 


=  868.5°fc 


T3-T2 


T5-T2 

T3    =    e{T5-T2)  +  T2 

T3     =     (.8)(868.5  -  486.8)  =  792.2°fc 


Similarly 


n    =    Ts-e{Ts-T2) 

Te    =    868.5  -  (.8)(868.5  -  486.8)  =  563.1°fc 

Since  this  is  a  counterflow  design  with  C*   >  0.5,  the  best  choice  for 
average  temperatures  is  the  arithmetic  average  on  each  side 

Ts+Te       868.5  +  563.1 


Thm      = 


=  715.8% 


J-  cm.       — 


T2  +  Tz       486.8  +  792.2 


=  639.5% 


2  2 

Inlet,  outlet  and  average  densities  are  next  computed.    Pv  =  RT  or 
V  =  RT/P  but  V  =  1//?,  so  p  =  P/RT 

1.126  X  lO^N/m^ 


PS 


{287j/kg-°  k){868.5°k) 


=  0A517kg/m' 
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1-126  X  10^  3 


1)  =     1  (1+1)=  1.82  WA^ 


5.63  X  10^         ^  ^^^,    ^    3 
^  287  (486.8  ^^ 


5.63  X  10^  ,     ;    3 

P3     = =  2.476A:g/m"* 

'^^  (287  (792.2)  ^' 


-J  =     O.Z260m^/kg 


»n,c 


B.2      A^,u 

For  a  couterflow  heat  exchanger, 

l-expl-Ntujl-c*)] 
l-c*exp[-Ntu{l-  c*)] 

but  for  the  balanced  flow  case  (c*  =  1)  under  design, 

Ntu 
'      l  +  Ntu 

or 

^-=137=1^^  =  ^ 

Since  gas  is  the  working  fluid  on  both  sides,  it  was  estimated  that  both  sides 
would  have  approximately  the  same  surface  resistance  so  that 

Nt^,h  =  Ntu,c  ^  2Ntu  =  8 

B.3      Estimated  Pressure  Drops 

Recall 

^"^  '''■'  =  1,1193 


J-  c  /  avg         boy. 5 
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and 


9  _  1.1193 
^~  (175)2" 


27.983 


APf,  _   k{9/r^)   _  (0.04) (27.983) 
Ph    ~  1  +  9/r^  ~      1  +  27.983 


0.03862 


AP. 


i  +  e/i 


0.04 
1  +  27.983 


=     0.00138 


From  Figure  2.3,  an  estimate  of  j/f  =  0.0035/0.0077  was  obtained. 


B.4     Mass  Velocities 


The  mass  velocity  on  each  side  is  estimated  from  the  relationship 


G  = 


\     Ntti     J  aide  f 


As  a  first  approximation,  t]o  was  assumed  to  be  0.90 

(2)(1.126x  10^ iV/m^) (0.9) \  /0.03862\  /0.0035\ 


Gh     = 


(1.824mVA:y)(0.7)3 


2 


8       /  V  0.0077/ 


Oh     =     I7.58ekg/m^  -  s 


G,    = 


(2)(5.63  X  105)(0.9) 
(.3260)  (0.7)  I 


0.00138\  / 0.0035 \ 
8       )  V  0.0077/ 


Gc     =     U.SSZkg/rn^  -  s 


NRe,h 
NRe,h 


GhDh       {ll.h^Zkg/m^  -  s)(0.0474/f)(.3048m//f) 


/^fc 


=     7481 


3.396  X  IQ-^kg/m-  s 
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=     (17.583)(.0474)(.3048)  ^ 
^*''  3.159  X  10-5 

Entering  Figure  2.3  with  these  values  yields 

fh     =     0.0078 

jh    =    0.00355 

fc    =     0.0077 

jc    =     0.0035 

B.5      Heat  Transfer  Coefficients  and  Fin  Effective- 
nesses 

With  the  Colburn  factor  known,  the  heat  transfer  coefficient  can  be  deter- 
mined from 

This  result  is  used  to  compute  the  fin  characteristic  length, 

The  temperature  eff'ectiveness  of  the  fin  is  next  calculated  from 

tanhml 


ni 


ml 

A  good  approximation  for  /  is  given  by  half  the  plate  spacing  minus  the  fin 
thickness. 

The  total  surface  temperature  effectiveness  can  then  be  calculated 

r/o  =  1  -  (1  -  r)f)Af/A 

where  A// A  is  a  characteristic  of  the  surface  selected. 

,  (0.00355)(17.586A:g/m2  -  s){1.0kj/kg  -°  k) 

hh     =     ^ 

(0.7): 


3 
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hh    =     79.19w/m^  -°  k 


rrih    = 


{2){79.19w/m^  -°  k) 


(190«;/m  -°  A:)(0.032m)(0.0254m/m) 


rrih    =    32.024m 


-1 


/     = 


0.750 


in  -  .032xn)  0.0254m/m 


/     =     8.71  X  IQ-^m 


^I,h     = 


fan/i(32.024)(0.71  x  IQ-^) 
(32.024)(0.71  X  10-3) 


Vf,h 


=    0.9749 


r,,,h    =     1- (1-0.9749)  (.606) 


h,     = 
rric     = 


=     0.9848 


(0.00355)(17.583)(1.0) 
(07)1 


(2) (78.06) 


=  78.06u;/m^  -°  k 


(190)(0.032)(0.0254). 


=  31.79m 


-1 


tan/i(31.79)(8.71  x  lO'^) 
(31.79)(8.71  X  10-3) 


=  0.9752 


r^o.c     =     1  -  (1  -  .9752)(.606)  =  0.9850 
The  overall  heat  transfer  coefficient  can  now  be  determined  neglecting  foul- 
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ing  resistances  and  wall  thermal  resistances, 

1  1  1 

+ 


1  1 


+  z r; r  =  0.02583 


U  (0.9848)  (79.19)       (.9850)  (78.06) 

U    =  Z8.7Uw/m^  -°  k 

B.6      Dimensions 

The  heat  transfer  surface  area  is  computed  from 
NtuC 


A    = 


U 


A    =     (4)(3600fcVH(l-0fc.7yfcg  -°  fc)  ^  103  30^2 
38.714u;/m2  -°  k 

Mininaum  free  flow  area  is  given  by 

A„  =  [W/G]  =        ^^^''3/h'        =  0.05686m2 

In  order  to  determine  frontal  area,  a  must  be  computed  from  the  relationship 

dfiDh/4       ■ 
d  +  d  +  2t 

Here  a  parting  plate  thickness  of  OAmm  was  assumed,  so 

(0.75»n)(76.1/fV/f3)(0.0474/4/0 


a     = 


0.75»n  +  0.75»n  +  2(.4mm)(. 1/2.54  in/mm) 
a    =     0.44162 

Afr      =      Ao/o- 


0.05686m2  , 

Afr    =     =  0.12876m2 

^"^  0.44162 
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The  flow  length  is  given  by 


L    = 


4Ac 

(0.0474/t)(0.3048m//t)(l03.321m2) 


(4)(0.05686m2) 

L    =     6.56m 
The  total  volume  of  the  heat  exchanger  is  calculated  by 
Vol     ={Afr)iL) 

m 

Vol  =  (0.12876m2)(6.56m)  =  0.845m^ 

B.7      Calculated  Pressure  Drops 

For  this  level  of  design,  entrance  and  exit  eflFects  are  not  considered,  so  only 
core  friction  effects  will  be  used. 

AP  .,,  ,    __  .     „,  ,_„ 


0.00125 


P 

j{ij/rh){i/p)m^  /^-r 

APh 
Pk 

= 

(0.0078)(6.56)(1.824)(17.586)2 

(2) (0.0474/4) (.3048) (5.63  x  10^) 

Similarly 

^  =  0.00125 

so  that 

AF 
S-— =  0.03674 

The  cycle  efficiency  with  these  input  parameters  is 

n 

(T4  -  Ts)  +  (Ti  -  T2) 
T4-T3 

n 

(1300  -  868.5)  +  (300  -  486.8) 

(1300  -  792.2) 
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ri     =     0.4819 
and  the  efficiency  degrade  is  given  by  7700  —  V- 

r}co-r}^-  0.5760  -  0.4819  =  0.0941 
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Appendix  C 

Minimum  Entropy  Design 
Method  Numerical  Example 


This  example  uses  the  data  from  case  C-12  of  Table  3.1, 
G  =  12A56kg/m^  -  s,  j  =  0.00355,  and  /  =  0.0078. 

C.l      Optimum  Length 

In  order  to  use  equation  (4.24),  intermediate  values  must  first  be  determined 

7  0.00355 

The  temperature  span  parameter,  r,  cannot  be  calculated  until  T5  is  known, 
so  an  estimate  based  on  previous  iterations  is  used.  If  the  calculated  value 
of  T  is  close  to  the  value  used,  the  computations  are  valid. 


{L/Th)opt      = 

2   /               r/Nst               Y 

G  \Bf{l/p,P,  +  l/pHPh)) 

{L/rh)opt    = 

2       f                     0.365/0.00450 

12.456  ^(o.287)(0.0078)  (,^M2|_,  +  _yi8_^) 

{L/rh)opt     = 

5283 
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Since  rh.  is  set  by  the  geometry  selected,  L  can  be  calculated 
Lopt    =    5283(.0474/f)(.3048m//0 

Lopt    =     19.08m 

C.2      Heat  Exchanger  Effectiveness  and  Pressure 
Drops 

From  equation  (4.16), 

Ntu,h    =     (5283)  (0.00450)  =  23.79 
Applying  the  same  assumption  as  was  used  in  Chapter  3, 

Ntu,h     =     Ntu,c  =  2^■tu 

Ntu    =    23.79/2  =  11.89 

Then,  since  this  is  still  a  balanced  counterflow  design, 

^tu  11.89  ^^„„ 

€  =  -i—  =  =  0.922 

l  +  Ntu       1  +  11.89 

For  hot  and  cold  side  pressure  drops,  equation  (4.17)  is  used. 

AP 
{~p-)H     =     f{L/n)g'/{2pP)  =  fiL/r,)g' 

AP 
(-—)/,     =     (0.0078)(5283)(0.05005)2  =  0.010322 

(-—)<:     =     (0.0077)(5283)(0.00946)2  =  0.00369 

S— -     =     0.10322  +  0.00369=0.1069 
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C.3      Entropy  Generation 

Equations  (4.14)  and  (4.15)  will  be  used  to  determine  entropy  generation  on 
both  sides. 


"''  =  Nz:^\7JA-F) 


0.365 
NsH    =     ^^^  +  (0.287)(0.10322) 

Nsh  =  0.04497 

Ns,  =  ^  +  (0.287)(0.00369) 

Nsc  =  0.0164 

Ns  =  Nsh  +  Nsc 

Ns  =  0.04497  +  0.0164  =  0.0614 

C.4      Heat  Exchanger  Volume 

Since  Afr  is  calculated  from  minimum  free  flow  area  and  other  surface  geo- 
metrical constraints, Ay^  from  the  corresponding  case  in  Chapter  3  was  used 
to  calculate  the  heat  exchanger  volume. 

Vol  =  L{Afr)  =  19.08m(0.12876m^)  =  2.46m^ 

C.5      Cycle  Efficiency  Degrade 

Temperatures  that  reflect  the  actual  effectiveness  and  pressure  drop  must 
first  be  determined. 

-^    =     -^(1  -k)  =  5(1  -  0.1069)  =  4.465 
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p  \'7pt("v"^ 


Ts    = 


ik) 


-\  -1 


Ts    = 


(4.465) 


(.9)(1.4-1) 


1.4 


1300 


=  884.8°A: 


Ts    =    e{Ts-T2)  +  T2 


Ts    =    0.922(884.8  -  486.8)  +  486.8 


Ta     =     853.8°A: 


V     = 


(T4  -  T5)  +  (Ti  -  T2) 
T4-T3 


(1300  -  884.8)  +   300  -  486.8 

T7     =     ^^ ^-—^ ^  =  0.5119 

'  1300  -  853.8 


r]^-rj     =     0.5760-0.5119  =  0.06412 
A  check  on  the  vaUdity  of  the  r  value  assumed  shows 

(  Ph     Fn\' 


T      = 


884.8         /486.8 


486.8       V  884.8 


r    =     0.367 


Therefore,  assumed  value  of  0.365  was  valid 
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Appendix  D 

Sample  Calculation  of 
Entropy  Generation  for 
Conventional  Designs 


Example  uses  case  C-12.  From  previous  example, 

T    =     0.34461 

4r^     =     0.0359 
Ph 

—A    =    0.00136 


Ntu,h    =    Ntu,c  =  8 
Using  equations  (4.14)  and  (4.15), 

T       [ r\   zap 

Ns.h      =       TT +       — 


Similarly, 


Ntu,h       \^pJh^  ^  ^  ^ 


Ns,h     = +  (0.287)(0.0359)  =  0.05338 

8 


Ns,c    =     ^-z +  (0.287)(0.00136)  =  0.04347 
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N,     =     Ns,h  +  Ns,c  =  0.09685 

Adding  the  first  term  of  each  side's  expression  gives  the  AT  contribution, 
and  the  sum  of  the  second  terms  is  the  AP  contribution  to  the  entropy 
generation. 

Ns,AT     =    0.08615 
N,  AP     =     0.01070 
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